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Abstract 

The technologically important problem of predicting Auger recombination lifetimes in semicon- 
ductors is addressed by means of a fully first-principles formalism. The calculations employ highly 
precise energy bands and wave functions provided by the full-potential linearized augmented plane 
wave (FLAPW) code based on the screened exchange local density approximation. The minority 
carrier Auger lifetime is determined by two closely related approaches: i) a direct evaluation of 
the Auger rates within Fermi's Golden Rule, and ii) an indirect evaluation, based on a detailed 
balance formulation combining Auger recombination and its inverse process, impact ionization, in 
a unified framework. Calculated carrier lifetimes determined with the direct and indirect methods 
show excellent consistency i) between them for n-doped GaAs and ii) with measured values for 
GaAs and InGaAs. This demonstrates the validity and accuracy of the computational formalism 
for the Auger lifetime and indicates a new sensitive tool for possible use in materials performance 
optimization. 

PACS numbers: 79.20.Fv,79.20.Ap,71.15.Mb 
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Led by advances in computational simulation, the paradigm within semiconductor de- 
vice engineering is shifting away from the static perspective of material properties as design 
"constraints" towards a more empowering view of material properties as design optimization 
parameters. This growing trend is evident clearly in recent initiatives in nano-engineered 
materials from quantum dots, quantum wells and nanowires to band structure engineering 
of conventional |I| and ordered superlattice materials. However, despite the fundamental 
importance of these processes, owing to a lack of adequate capability for their prediction, 
carrier recombination processes in semiconductors have remained largely beyond the reach 
of materials design. Minority carrier lifetime is a critical, performance-limiting material 
parameter in many opto-electronic devices, as well as in bipolar transistors, p — n rectifiers 
and thyristors. Two recombination processes - radiative and Auger - generally impose the- 
oretical and practical limits on carrier lifetimes. Minority carrier lifetimes in heavily doped 
and/or narrow band gap materials (such as 1.6 fim Ino.53Gao.47 As - a basic material for 
thermophotovoltaic devices) tend to be limited by Auger recombination (AR), and so its 
first -principles determination is the focus of this paper. 

In n-doped, direct band gap materials, the dominant Auger process is electron-electron- 
hole (e^e^h) recombination [0. As shown in Fig. ?? (a), in the e^e^h process, a valence band 
hole decays via non-radiative recombination with a conduction band electron, conserving 
energy and momentum through transfer to a second conduction band electron. This process 
commonly is referred to as the conduction-conduction-conduction- "heavy hole" (CCCH) 
process. In the inverse process of electron-initiated impact ionization (I-I), depicted in Fig. 
?? (b), an energetic primary electron promotes a secondary valence band electron into the 
conduction band, producing a mobile electron-hole pair. Processes closely analogous to 
those depicted in Fig. ?? (a) and (b) occur in p-type materials, resulting in an electron 
decay through a "conduction electron" -"heavy hole" -"heavy hole" -"light hole" (CHHL) 
AR and in a hole-initiated I-I, respectively. 

In this work, the fully first-principles approach that was used successfully to evaluate I-I 
rates is extended to the calculation of AR rates which are found to be in good agreement 
with previous experimental and theoretical values reported for the most studied semicon- 
ductor - n-type GaAs - and for technologically important InGaAs. Moreover, AR lifetimes 
calculated directly via Fermi's Golden Rule are shown to be strikingly consistent with life- 
time values determined via detailed balance principles proceeding from calculated I-I rates. 
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To our knowledge, this density functional theory (DFT[0) study represents the first fully 
ab initio determination of Auger recombination lifetimes. The agreement of computational 
results with experiment for n-type GaAs and Ino.5Gao.5As indicates that this method is 
sufficiently precise and reliable to be able to correctly predict trends of the Auger coefficient 
as a function of physical parameters, such as doping, composition, etc. Since the current 
approach is all-electron and ab initio in nature, we expect the same methods to be equally 
successful for new and more complex materials. Finally, we find that detailed balance eval- 
uations proceeding from I-I rate calculations appear to be a computationally efficient means 
of determining AR rates. 

Within an independent particle scheme 0, given a hole in state (n2,k2), the AR rate is 
expressed as: 

fiKMKimElll + K: - Kl - Kl) (1) 

and the total AR rate is: 

^^"^ = E / g{E^) R^^'in,, k^) (2) 

Here, the rii are the band indices {ni,i = 1,...,7 {i = 5,.. .,11) for GaAs (InGaAs), corre- 
sponding to the four valence and three conduction bands considered) and ki 3 4 are k points 
in the full Brillouin zone (BZ) |^. The many-fold integration over the BZ is performed 
according to the Sano-Yoshii scheme]^, using 152 and 126 k points in the irreducible 
wedge of the zincblende and tetragonal Brillouin zone, respectively [l^. The summation 
over k2 (Eq.^ is carried out in the irreducible wedge of the zone. Energy conservation is 
enforced through the S function in the energy eigenvalues; / denotes the Fermi-Dirac oc- 
cupation probability and g = 1 — f. The matrix elements, taking into account both direct 
and exchange terms (obtained from the direct contribution through the exchange of final 
states), are of the form: \M\^ = |(|MdP + \Me\'^ + \Md - Me\'^). The direct term is given 
as: 

Md = ^ ^(ks + k4 - kl - k2 + Go) (3) 
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Pn3,k3;ni,ki n4,k4;n2,k2 (Go — Gu) 

£(g)(|ki-k3 + GuP + A2) 

where e is the electronic charge and Q is the volume of the unit cell. Momentum conservation 
is enforced through the S function in the k vectors and Pn^ ,kf;ni,ki(G) is the Fourier transform 
of the overlap matrix of the wave functions, i.e., p„^,kf;ni,ki(r) = ^n^,kf ('^)^ni,ki(i")- The 
subscripts i and / denote initial and final states; q = |ki — ks + Gu| is the momentum 
transfer and Go and Gu are reciprocal lattice vectors. To describe the interaction between 
the valence and conduction band electrons, we used a Coulomb interaction screened through 
a static dielectric function proposed by Cappellini et ai, 6{q) ||lT|, that has been shown to 
be very accurate for semiconductors. The interaction between conduction band carriers is 
modeled through a Debye potential (screened Coulomb interaction), in which the inverse of 



the screening length is expressed as A = y^^^, T and Kb being the carrier density, 
the temperature (here, T=300 K) and the Boltzmann constant, respectively. 



The theoretical estimate of the I-I rate|]12|, |13| is known to be quite sensitive to the 
assumed band structure. Therefore, the importance of employing accurate quasi-particle 
bands as a basis for I-I and AR calculations needs to be emphasized. To faithfully re- 
produce experimental band structures (both occupied and virtual states), we employ the 
screened-exchange local density approximation (sX-LDA) |T5[, as implemented ||1 6| | self- 
consistently in the highly precise, all-electron, full-potential linearized augmented plane 
wave (FLAPW) method. The sX-LDA approach corrects known deficiencies of the LDA 
for excited states and has been shown to be very successful for treating a large variety of 
semiconductors [TgH . 



The effect of spin-orbit coupling, neglected here, was investigated in Ref. |T3[ for GaAs; 
the overall effect was shown to be relatively small. Umklapp processes for Gu = were 
fully included in our calculations; the Go 7^ terms, which are important at high energies 



were not included, since Auger is a threshold process. Phonon-assisted transitions [0, |T9 



and many-body|20, 21| effects were also not included; for InGaAs, Dutta and Nelson |19 



showed that phonon assisted AR should not be relevant compared to the direct AR. Within 
these approximations, giving an estimated overall error of 30-40 % on the AR rates, we show 
below that our simple band-to-band AR within DFT correctly describes the physical trends 
as a function of several technological parameters, and so may be used as a sensitive tool in 
materials device optimization. 
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It can be shown for band gaps Egap » KbT, that the Auger partial hfetime can 
be expressed as 



- = — (4) 

Ip ^'Ip 

where is the number of minority carriers per unit cell at equilibrium {Np/Q corresponds 
to the minority carrier density p)P^. In Fig.?? we show the dependence of the calculated 



hole Auger lifetime on carrier density, which is approximately linear on a log-log scale. 
Similar trends have been reported, for example, in an empirical pseudopotential study of 
AR in Si and in a recent experimental study of AR in n-doped InGaAs p3| . 
The Auger coefficient, C„, is defined by the relation: 

- = C^iny (5) 

Recall that partial recombination lifetimes usually are "determined" experimentally based 
on curve fits to carrier lifetime measurements resulting from several competing recombina- 
tion processes. The parabolic relationship above is the usual basis for determinations of the 
Auger recombination coefficient. It is noted, however, [^] that both shifts in the Fermi en- 
ergy (particularly pronounced in many III-V alloys with low effective masses), and changes 
in the Debye screening length introduce a carrier density dependence into Eqs.[^ and H, so 
that Cn itself is expected to depend on n. This is actually the case from our calculations. 
Indeed, in GaAs our C„, determined as the ratio of the inverse lifetime and the carrier 
density squared, ranges from 2.3x10^'^'^ cm^^s^^ at n ~ 4x10^^ cm^'^ to 0.5x10^'^° cm^^s^^ 
at n ~ 6x10^^ cm~^. This range of values is compared in Table |I| with previously reported 
experimental p^, |25| , |26[] and theoretical p7|, |28| results. In InGaAs, there is only a slight 



dependence of C„ on n, Cn ranging from 0.8x10 cm ^s ^ to 1.1x10 cm ^s ^ for con- 
centrations ranging from 7x10^^ cm~^ to 2x10^^ cm~'^. These values are compared in Table 
HI with previous experimental ||2^, ^ and theoretical [^, ^ results. The comparison 



between the calculated values for the binary and the ternary compound shows that there is 
a strong reduction (by an order of magnitude) of C„, suggesting that the Auger process is 
crucial in the narrow-gap InGaAs. Despite the wide range over which the reported Auger 
coefficients vary, the calculated trend with Ga content is confirmed experimentally, unam- 
biguously emphasizing the power of DFT in predicting both quantitatively and qualitatively 
correct physical trends as a function of doping, pressure, composition, etc. The wide range 



of reported Auger coefficients suggests the complexity of both the measurements and the 
calculations. Our result for the Auger coefficient is in good agreement with most of the pub- 
lished experimental values and is again consistent, when compared, with previously reported 
results from model calculations incorporating various simplifications (such as constant effec- 
tive masses [ P7| , All considered, we conclude that our results confirm the reliability of 
the computational methods and band structures herein described i) to accurately calculate 
Auger recombination lifetimes that cannot be easily or unambiguously determined by ex- 
periment; ii) for mapping trends with respect to material composition and dopant density, 
that are useful to guide optoelectronic device and materials design; and Hi) for describing 
Auger recombination in novel and more complex materials. 

Auger recombination and impact ionization are related as inverse microscopic processes 
through the principle of detailed balance (PDB), in the same manner as optical emission 
and absorption. In fact, in the radiative (non-radiative) case, the e — h generation process 
is the absorption (I-I), in which the initial object is the photon (highly energetic impacting 
electron); correspondingly, their inverse recombination process is emission (AR), in which 
the final product is the photon (Auger electron). 

Following Landsberg0, the AR lifetime is given by, 

- ~ 1^ E / d'k^P{n,, k,)R'-\n,, ki)/(ni, ki) (6) 

where /3(ni,ki) is a factor (< 1) that takes into account the occurrence of processes from 
state {ill, ki) other than I-I and the R^'^{ni, ki) are the I-I rates. This last expression is the 
non-radiative equivalent of the van Roosbroeck-Shockley relation |3^, that has been used 



extensively in determinations of radiative recombination lifetimes from the complex dielectric 
function. However, while experimental dielectric function data (or equivalently, refractive 
index and extinction coefficient data) are abundant and highly reliable, experimental data 
for electron-initiated impact ionization rates are normally inadequate to estimate the AR 
lifetime from Eq. (|^). A utility of the current formalism is that one can calculate both 
I-I rates and AR rates either independently from first principles, or from one another via 
detailed balance arguments. In our formalism, we focus on two selected processes (I-I and 
AR), each being the inverse of the other; therefore, the PDB holds with /?(ni,ki) = 1 (or, 
equivalently, that every energetic electron in state (ni,ki) initiates an impact ionization 
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process). In Fig. ??, Eq. (H) is applied to the previously calculated I-I rates (see 0] for 
details) and the Auger lifetime for = 1.0x10^^ cm~'^ is derived. The open symbol in Fig. 
?? denotes the value obtained, which is in striking agreement with the values calculated 
using the direct approach (see Eq.|ip. The excellent agreement achieved between the direct 
and "indirect" calculations of both parameters provides confidence in the equivalency of the 
two approaches and in the numerical procedures embodied therein. 

In summary, we have presented a fully first-principles formalism for calculating Auger 
recombination lifetimes using screened exchange FLAPW quasiparticle wavefunctions and 
band structures determined self-consistently. The numerical accuracy of the approach was 
checked using the principle of detailed balance applied to the two simulated processes, I-I 
and AR. The Auger lifetime was determined using two equivalent approaches, "direct" and 
"indirect," with highly consistent results. Our results for GaAs and InGaAs are in excellent 
agreement with the most accurate experimental and theoretical data, providing confidence 
in the computational method and justifying its future applications to more complex systems. 

Work at Northwestern University supported by the National Science Foundation (through 
its MRSEC program at the Materials Research Center). 
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TABLE II: Auger coefficients in n-type InGaAs. 



Exp. technique or theor. approach 


Cn (cm^s-i) 


sX-LDA FLAPW 


(0.8 <Cn< l.l)-10-29 


Time-resolved photoluminescence decay" 
Radio-frequency photoconductive decay^ 
Photoluminescence photon-counting 
+ degenerate carrier conditions 


0.5-10-29 

7-10-29 
18-10-29 


non-parab. bands + phonon-assisted + efF. masses'^ 
Kane model for wave functions'^ 


0.5-10-29 

3-10-29 
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h. Ref. 30 



c. Ref. 23 



d. Ref. 31 



e. Ref. 32 



Simulated processes in n-type materials: (a) Auger recombination, resulting in a hole- 
decay; (b) electron initiated impact ionization, resulting in a pair production. 

Calculated hole Auger lifetime vs carrier concentration in (a) GaAs and (b) InGaAs. The 
results for the binary compound are obtained according to the direct approach (see Eqs.(2-5) 
- filled circles) and to the indirect approach (see Eq.(6) - open circle). 
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